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Abstract 

A definition of K41 scaling law for suitable families of measures is given and 
investigated. First, a number of necessary conditions are proved. They imply 
the absence of scaling laws for 2D stochastic Navier-Stokes equations and for the 
stochastic Stokes (linear) problem in any dimension, while they imply a lower 
bound on the mean vortex stretching in 3D. Second, for 3D stochastic Navier- 
Stokes equations necessary and sufficient conditions for K41 are proved, translat- 
ing the problem into bounds for energy and enstrophy of high and low modes 
respectively. The validity of such conditions in 3D remains open. Finally, a 
stochastic vortex model with such properties is presented. 

1 Introduction 

In very rough terms, the scaling law devised by Kolmogorov and Obukhov for 
turbulent 3D fluids (usually referred as K41, see llT2l and a detailed discussion 
in [9 1), says that S^r) ~ e 2 / 3 r 2 / 3 where S^r) is the second order structure func- 
tion and e is the mean energy dissipation rate. Moreover, it is specified that this 
law is valid at very high Reynolds numbers and for distances r in a certain range 
between the integral range and the Kolmogorov dissipation scale, having the order 
rj = i/ 3 / 4 e -1 / 4 . Although the numerical essence of these claims may be clear, their 
precise mathematical interpretation is not necessarily unique and could change a 
little bit depending on new discoveries. 

The purpose of this note is to give one possible precise mathematical formula- 
tion of this scaling law and to discuss it from a number of viewpoints. We imme- 
diately stress that we cannot prove its validity for the 3D Navier-Stokes equations, 
but nevertheless we obtain a number of insights that seem worth to be known. 

Some of our considerations are true for quite general families of probability 
measures; others will be specific to the stochastic Navier-Stokes equations on the 
torus [0, l] d , d = 2,3, 

du 

-^ + {u-V)u + Vp = uAu + h a (x)/3 a (t) (1.1) 

a 
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with divu = and periodic boundary conditions, with suitable vector fields h a {x) 
and independent Brownian motions (3 a (t) (the torus instead of a more realistic 
framework has been chosen for mathematical simplicity). We consider this equa- 
tion in the limit v — ► 0. Since the force does not vanish as v — > 0, this is a singular 
limit problem much like the boundary layer one, and so may be considered as a 
prototype of high Reynold number singular limit problem, with some mathemati- 
cal simplification due to the advantages produced by stochastic analysis. It should 
be noted that another possible and interesting approach to the zero-viscosity limit 
is the one adopted in Fl3ll (for the 2d case), where the amplitude of the forcing 
noise is proportional to the square-root of the viscosity. 

In the following we shall use the parameter in place of the Reynolds num- 
ber; this simplification is justified in our model since the force and the domain are 
given, so the Reynolds number goes to infinity if and only if v — > 0. 

We shall denote by H the natural space of finite energy velocity fields on the 
torus and we shall introduce a space V of probability measures on H having certain 
symmetries and regularities (precise definitions are given in the next section). On 
the fields <p( l \x) we shall assume conditions such that there exists at least one 
stationary probability measure p, G V associated to (11.11) (stationary measures will 
be defined in the next section). We use the notation 



for some coordinate unitary vector e, with r > (the results proved below extend 
to the so called longitudinal structure function; we consider dl.2t to fix the ideas). 
The measures of V are supported on continuous vector fields, so the point wise 
operations in (11.21) are meaningful. Moreover, the symmetries in V imply that 
S*2 (r) is independent of the coordinate unitary vector e (in addition most of the 
estimates proved in the sequel extend to every unitary vector e). 

We are going to define K41 scaling law for a set M. C V x M.+. The rea- 
son is that equation (11.11) may have (a priori) more than one stationary measure 
for any given v and in certain claims it seems easier to consider a set of mea- 
sures for a given v. Given v > we use the notation M. v for the set section 

Given (p,u) &V x K + , we define the mean energy dissipation rate as 




whenever the integral is well defined. 

For every p G V we introduce the second order structure function 



S£(r) = M ||n(r-e)-u(0)|| 2 



(1.2) 



e = e(/i, v) := u ■ fx 



j 



[0,1] 



Du{x)\\ 2 dx . 



Remark 1.1 If fi is a stationary measure of ( li.7l) and a mean energy equality 
(coming from ltd formula) can be rigorously proved, one can show that e does 
not depend on (/x, v). 
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Given (//, v) € V x R+, we also define the quantity 

rj = Ti(ji,v) := v 3/4 e(fi,v)~ 1/4 . 

Remark 1.2 case 0/ equations ALU . r\ is a length scale: v has dimension 
[L] 2 [T] _1 , e has dimension [L] 2 [Tp 3 , 50 r\ has dimension [L\. The only com- 
bination of v and e in powers, having dimension [L], is the r/ above. This is the 
simplest reason to choose rj as a length scale involved in K41 theory. More refined 
arguments may be found in and related references. 

Let us come to the definition of K41 scaling law chosen in this work. Here and 
in the sequel, when we talk about a set M C V x IR + , we tacitly assume that 

M. v 7^ (j) for all sufficiently small v > , 

since otherwise several definitions and statements would be just empty. 

Definition 1.3 We say that a scaling law of K41 type holds true for a set Ai C 
V x M + if there exist vq > 0, C > c > 0, Co > 0, and a monotone function 
Rq : (0, uq] — > 1R+ with Ro(v) > Cq and lim^o Roi^) = +00, such that the 
bound 

c-r 2 ^ < S%(r) <C-r 2 / 3 
holds for every pair (fi, v) G M and every r such that v £ (0, vq\ and 

C ■ 7](fi, v) < r < rj(fi, v) ■ Rq{v). 

Remark 1.4 For simplicity we could have asked the scaling property for Cq ■ 
ijdijv) < r < rofor a constant tq (a measure of the integral scale). However, 
such a formulation could be too restrictive. On the other hand, it is necessary that 
the range ofr's increases to infinity (relative to ij) as v — > 0, otherwise the property 
becomes trivial, see remark \T3\ 

This is the mathematical formulation of K41 theory that we analyse in this note. 
Here is a list of facts we can prove around it. In summary, they have the structure 
of certain necessary conditions for K41, and certain almost equivalent conditions. 

• We introduce a measure 6 of the length scale where dissipation takes place, 
defined as 

^ J[o,i] d 11-°^) 



(1-3) 



^[/[o,i] d W D2u ( x )\\ dx 
and prove a claim of the following form (theorem l2.2t : if a scaling law holds 
on a range Co ■ r] < r < i] ■ Ro(v), then 

6 < Cr] . 



In other words, the scale at which dissipation dominates cannot overlap with 
the range over which a fractal scaling law holds. 
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Since 9 is constant (with respect to u) for both the 2D stochastic Navier- 
Stokes equations and the Stokes (linear) equations, we can rule out K41 
scaling law for such systems. For the theory of 2D stochastic Navier-Stokes 
equations this seems a remarkable fact. Moreover, these facts tell us that, in 
case definition ! 1 . 3l holds true in the 3D case, it is strictly due to 3D nonlinear 
effects. This is further emphasised by the following result. 

For the 3D stochastic Navier-Stokes equations we prove that, if K41 holds, 
then the mean vortex stretching 



| (S u curl u, curl u)\ 2 dx 



(where we set S u = ^(Du + Du T )) must be very large, essentially at least 
as large as z/~ 3 / 2 . See Corollaries I2.19l and l2.21l Vortex stretching is thus a 
basic mechanism in K41 theory. 

• We apply a known scaling transformation (see (0) and introduce an auxiliary 
family l3~2l of stochastic Navier-Stokes equations with modified domain and 
viscosity. Then we introduce a condition on this family of equations, called 
Condition^ and prove it is equivalent to the scaling law of K41 type. This 
is conceptually interesting since Condition |A] is of rather qualitative nature, 
while its consequence (the scaling law of K41 type) is more quantitative. 
More specifically, the behaviour r 2 / 3 , and also the exponent 3/4 in the def- 
inition of rj, arise from the stochastic Navier-Stokes equations themselves, 
through a scaling transformation, when certain bounds (without special ex- 
ponents) are fulfilled for the auxiliary family. 

• We give other necessary and sufficient conditions for K41, starting from 
Condition |A] In plain words, they state that the large structures of u have 
bounded mean square gradients (or bounded mean ens trophy), while the 
small structures have bounded mean energy. 

• As a mild support to the belief that all these Conditions could be true in di- 
mension d = 3, we finally exhibit a random field that satisfies them, and 
was constructed independently from this purpose in Q as a model of tur- 
bulent fluid inspired by the vortex structures usually observed in numerical 
simulations. 

The remainder of the article is structured as follows. The remainder of the 
present section introduces the notations that will be used throughout this work. In 
Section |2j we draw several conclusions from our formulation of the K41 scaling 
law that allow us to give stringent necessary conditions for it to hold. These con- 
ditions are sufficient to rule out any non-trivial scaling law in the 2D case. We 
proceed in Section|3]to find a condition that turns out to be equivalent to K41. This 
condition is then shown in Section|4]to hold for a random eddy model introduced 
inQ. 
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1.1 Notations about functions spaces 

Let T be the torus [0, l] d , d = 2, 3, L 2 (T) be the space of vector fields u : T -> R d 
with L 2 {T) -components, HP(T) be the analogous Sobolev spaces, C(T) be the 
analogous space of continuous fields. 

Let H be the space of all fields u G L 2 (T) suchthatdivu = and J r u(x)dx = 
(zero mean) and the trace of u ■ n on the boundary is periodic (where n is the outer 
normal, see fl31 . Ch.I, Thm 1.2). Let V be the space of divergence free, zero mean, 
periodic elements of H 1 (T) and D(A) be the space of divergence free, zero mean, 
periodic elements of H 2 (T). Finally, let V be the space of infinitely differentiable 
divergence free, zero mean, periodic fields on T. The spaces V, D(A) and V are 
dense and compactly embedded in H. Let A : D(A) C H — > H be the (Stokes) 
operator Au = — An (componentwise). 

Sometimes we shall also need the same framework for the torus [0,L] d , d = 
2,3, with any L > 0. We set T L = [0,L] d , H L equal to the set of all fields 
u G L 2 (7i) such that div u = and u • n on the boundary is periodic, Vl, D(A£) 
and A L : D(A L ) C H L — > the analogs of V, and A. Notice only that 

we define the inner product as 



1.2 The class V of probability measures 

Let Vo be the family of all probability measures \i on H (equipped with the Borel 
<T-algebra) such that fi(D(A)) = 1 (D(A) is a Borel set in H). Since H 2 (T) C 
C(T) by Sobolev embedding theorem, the elements of D(A) are continuous (have 
a continuous element in their equivalence class). Consequently, given xq G T, the 
mapping u ^ u(xq) is well defined on D(A), with values in M. d . In particular - , any 
expression of the form 

V[f(u(xi),...,u(x n ))} 

is well defined for given x\, ...,x n G T, given u G Vo, and suitable / : M nci — > 
E (for instance measurable non negative). It follows that S% (r) is well defined 
(possibly infinite) for every u G Po- 

The same argument does not apply to Du(xq) and D 2 u{xq), at least in d = 3. 
This is why we use lengthy expressions like 




(So that, roughly speaking, \u 



n(0)| 2 for homogeneous fields.) 



X 



Du(x)\\ 2 dx , u 



X 



-D 2 it(x)|| 



which are meaningful (possibly infinite) for every fi G Po- 
We denote by "P the class of all u G "Pq such that 
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and, for every a G T and every rotation R that transforms the set of coordinate 
axes in itself, 



fi[f(u(- - a))} = fi[f{u(R-))] = n[f(Ru(-))} 



(1.4) 



for all continuous bounded / : H — > R. In plain words, we impose space ho- 
mogeneity and a discrete form of isotropy (compatible with the symmetries of 
the torus). In the following we will refer to this symmetry as partial or discrete 
isotropy. 

Discrete isotropy is imposed for two reasons. On one hand, S% (r) is indepen- 
dent of the coordinate unitary vector e, since given two such vectors e, e' there is a 
rotation R as above such that Re' = e, so 



(i 



\u(re) -u(0)|| J 



\u 



{Rre) - u(R0)\ 



\R(u(re) -it(0))| 



|ii(re') — u(0)\ 



On the other hand, we use discrete isotropy through lemma lA3l in Appendix 1. 
Finally, notice that 5^ (r) < oo for every r > and p, G V, by Lemma !2.1l below. 



2 Necessary conditions for K41 
2.1 General results 

The results of this subsection apply to suitable families of probability measures, 
without any use of the Navier-Stokes equations. They will be applied to the stochas- 
tic Navier-Stokes equations in the next subsection. 

Given a measure \i € V, 7^ So, we introduce the number 9 = 9(fi) defined 

= 00. If fj, = 5o, 
1. We have 9<C 



by the identity (11.3b . letting 9 = when fi f T \\D 2 u(x)\\ dx 
numerator and denominator vanish and we arbitrarily define 9 = 
where the constant is universal and depends only on the Poincare constant of the 
torus. By the definitions, we have 



e{n,v) 



f T \\D 2 u(x)\\ 2 dx 



for every pair (fi, v) G V x R + . 

When, as in our application, the elements u G H have the meaning of velocity 
fields, by dimensional analysis we see that 9 has the dimension of a length. We 
interpret it as an estimate of the length scale where dissipation is more relevant. 
Indeed, very roughly, from 



j T \\D 2 u(x)\\ 2 dx 
f r \\Du(x)\\' 2 dx 



E|fc[ 2 (|fc| 2 |£(fc)l 2 
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we see that 9(fi) has the meaning of typical square wave length of dissipation 
(looking at |A;| 2 |{i(/i;)| 2 as a sort of distribution in wave space of the dissipation). 

Lemma 2.1 For every (j, G V such that 8(fi) > Owe have 



— ■r l < 
4d 



j r \\Du(x)\\ 2 dx 



(2.1) 



for every r E (0, 9{p)/Ad\. The upper bound is true for every r > even if 
d{p) = 0. 

Proof. Since we want to use Taylor formula for elements of D(A), we use the 
mollification described in Appendix 1. We denote by fi £ the mollifications of fi. 
We prove in Appendix 1 that, for given r and fi, 



limu e \\Du(0)\\' 

e— >0 L 



lim u £ 



|L> 2 n(0)| 



= A» 


[X 






= A» 


[X 







|n(re) — ii(0)|| 2 
By space homogeneity of \i e 



lim u e 



dx 

||l) 2 u(x) || 2 <ix 
|u(re) — ii(0)|| 2 



\u(re) — u 



Me 

\\Du 



Du(ae) 

2 



da 



and thus, by the previous convergence results, 



\u(re) - u(0)\y 



< r 2 n 



\\Du(x)\\ 2 dx 



T 



This implies the right-hand inequality of d2.lt for every r > 0. 
On the other hand, for smooth vector fields we have 

u(re) — u(0) = Du(0)re + r 2 D 2 u(ae)(e, e)da 



and thus 



Me 



\Du(0)re\ 



u(re) — it(0)||' 
l 



D u(ae)(e, e)da 
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Again from space homogeneity of p, £ , 



r 2 / D 2 u(ae)(e, e)da 



and from lemma lA~3l of Appendix 1 

ji e \\Du{ti)e\\ 2 

Therefore 



2" 








, 4 

<r fi £ 


~\\D 2 u(0)\f 



\Du{Q)\\' 



\u(re) — it(0)||" 



r 

> — ;Me 
- 2cT 



\Du(0W 



4 

r He 



|-D 2 u(°)| 



We thus have in the limit 



T 



| Du(x) || cfa; 
and therefore, by definition of 



4 



|-D 2 u(x)|| dx 



■T 2 . 



\\Du(x)\\ 2 dx 

This implies the left-hand inequality of d2.lt for r G (0, ^r]. The proof is com- 
plete. ■ 

Theorem 2.2 Le? C f x IR + &e a set with the following scaling property: 
there is a function rj : Ai — > R+ (the length scale of the scaling property), a 
decreasing function Rq : [0, oo) — > W.+, with lim^o Rq{v) = +oo, a scaling 
exponent a G (0, 2) and constants C2 > C\ > 0, C3 > 0, v$ > 0, such that 
Ro(v) > C3 and 



C x -r a < S%{r) <C 2 -r a for r G [C*n{n, v)Ro{y)] 



(2.2) 



for every v G (0, Uq) and every \i G M u . Let 0(p) be the dissipation length scale 
defined above. 

Then the two length scales 9{p) and rj(/i, v) are related by the property 



lim sup sup 



( 



9(p) 







< 00. 



(2.3) 



Proof. It is intuitively rather clear that d2.lt is in contradiction with (12.2b if 
the ranges of r where the two properties hold overlap, so we need the bound (12.3b - 
The proof below confirm this intuition by ruling out the possibility that the factor 

p, J r \\Du(x)\\ 2 dx may produce a compensation. 
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We argue by contradiction and assume that there exists a sequence (fi n , v n ) G 
M., with v n — > 0, such that 

lim Q^n) = +oo> (24) 

Notice that, in such a case, 6(fi n ) must be positive, so lemma |2~T1 applies. Let us 
consider two sequences r' n and r" defined as follows: 

r n = C 3 ??(^«> *V>, ^ = r' n a n 

with 

lim a n = +00, r" < r}(/j, n , u n )R (u n ), r" < ^ 



4d 

where we ask that the last two inequalities are satisfied at least eventually. Such a 
sequence r" exists because lim^o Ro(v) = +00 and i2Al is assumed. 

We have (eventually) r' n , r£ G (0, and r^, G [C 3 rj(fX n , u n ),r}(fj, n , v n )Ro{^ 

hence for both r n := r' n and r n := we have 

CiC < ^" {rn) < C 2 <, < (r n ) < f3 n r 2 n 

where we have set f3 n = fi n j T \\Du{x) \\ 2 dx . The contradiction will come from 

the fact that, if it could happen that (3 n adjusts the factor r\ to produce r", this 
cannot happen simultaneously for the two sequences r n = r' n and r n = r l ! a . Indeed, 
from the previous inequalities we must have 

C ir %<[3 n rl {3 n r 2 n <AdC 2 r% 

hence 

for both r n = r' n and r n = r". But the inequalities 

f3 n >C 1 {r' n ) a ~ 2 , /5„<4dC 2 «)^ 2 
and the assumption a < 2 imply 

r' n > Crl 

eventually, for a suitable constant C > 0. This is impossible since lim n _>oo a n = 
+00. The proof is complete. ■ 



Remark 2.3 The divergent factor R$(v) in the definition \2.2l of a scaling law is 
essential to have a non trivial definition. If on the contrary, we simply ask that the 
scaling law holds on a bounded interval r G [C^i],^, C^r] v \, we have a definition 
without real interest. Let us explain this fact with a (useless) definition and an 
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example. Let us say that a family A4 C V x R + satisfies a local a property, a < 2, 
if there is a function rj(p, v) and constants C2 > C\ > 0, C4 > C3 > 0, Vq > 0, 
such that 

C ir a < Sg(r) < C 2 r a for r G [C 3 v(fi, v),Ctf{p., v)} (2.5) 

for every v G (0, z^o) every \x G Al^. As an example, consider a case with the 
mapping v 1 — > M. v which is single valued and injective and 

Sf(r) = u- l r 2 

where Ai u = {fJ- u }- This function S% (r) certainly does not have any interest- 
ing scaling exponent (different from 2) but satisfies the previous local a property 
simultaneously for a continuum of values of a. Indeed, given any a G (0, 2) 
take rj(fi u ,u) = v 2 ~ a ; then given a choice of C4 > C3 > 0, for every r G 
[C3?j([j 1 u , v),C4rj(fi u , u)\, namely for v'^^r G [C3, C4], we have 

Sf(r) = (is-^ r y~ a r a G [Ci,C 2 ] -r a 

with C\ = C| _a , C2 = Cl~ a . This example shows that the local a property is 
not a distinguished scaling property. Moreover, it shows that i2.lt and ( 12. 51 1 are 
compatible: this is why a proof of theorem \2. 21 is necessary. 

Example 2.4 Since we have just given a negative example ( artificial, but close to 
what happens in 2D), let us also give an example of a function of '{y, r) which sat- 
isfies the properties of deftnition \1.3\ and also X2.1\ (to see that they are compatible). 
It may look artificial, but it was devised on the basis of the vortex model of $7$, 
described also below. The function is 

with rj = i^ 3 / 4 . We have 

r<^5f(r)= /V 3 (y) 2 f = \r*[u^-l] 
which is essentially the behaviour 12. 71 On the other hand, 

-mi^>>=/>H>(t)t 

= ^ r 2 /3 _ ^ 1/2 _ 3 2 

4 2 4 

which is bounded above and below by the order r 2 / 3 since r G [V 3 / 4 , l] (V 1 / 2 < 

r 2/3| 
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Let us finally state two general consequences of the previous theorem, that we 
shall apply to stochastic Navier-Stokes equations. 

Corollary 2.5 Given a family M C V x R + , if 



inf 0(p) > 

then no scaling law in the sense of the previous theorem may hold true with a length 
scale rj{fi, u) such that 



liminfj inf r](n,u) ) = 0. 



We shall see that this simple corollary applies to the 2D stochastic Navier- 
Stokes equation and the Stokes problem, so K41 scaling law is ruled out for these 
systems. 

Let us apply the theorem to the case of K41 scaling law. We take, in the previ- 
ous theorem, 

rj([i, v) = rj(fi, v) = u 3/4 e(fi, vy 1,A 



as in the introduction. In the following result, (i 
nite. 



J r \\D 2 u(x)\\ 2 dx 



may be infi- 



Corollary 2.6 Let M cPx M + be a family with the K41 scaling law, in the sense 
ofDeftnition \1.3\ Then there exist vq > and C > such that 



J \\D 2 u(x)\\ 2 dx >Ce 3 / 2 {n,u)-u~ 5 / 2 

for every v £ (0, vq) and every \i E M. v . 

Proof. From (I2.3t . the definition of r)(fj,, v) and the definition of 9 2 (fi) we have 

limsup sup — — — F : — r | < oo. 





j T \\Du(x)\\ 2 dx 




u s ^ 2 e(fi,u) X ^ 2 /i 


J r \\D 2 u(x)fdx 



Thus, from the definition of e(/x, v), 



lim sup sup 



v 5 / 2 e(/i,i/) 



3/2 



< OO. 



v-+o~ \ueM„ p J r \\D 2 u(x)\\ 2 dx 
This implies the claim of the Corollary. ■ 

Remark 2.7 Dimensional analysis says that v has dimension [L] 2 [T]~ X , e has 
dimension [L] 2 [T] -3 , so e^/ 2 {^,u) ■ i/~ 5 / 2 has dimension [L]~ 2 [T]~ 2 , the correct 

dimension of J T \\D 2 u(x)\\ 2 dx 
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2.2 Application to stochastic Navier-Stokes equations 

In this section we consider equation (ll.lt in dimension 2 and 3 and also the corre- 
sponding linear equations (Stokes equations). 

2.2.1 The noise 

Since we are dealing with spaces of translation invariant measures, we wish to 
consider classes of noises that produce such measures. Every Gaussian translation 
invariant noise is 'diagonal' with respect to the Stokes operator A in the sense 
that eigenmodes are all independent. In order to give a rigorous definition for our 
driving noise, we define 

A (oo) := h E 2TrZ d : \k\ > Cjj 
and we assume that the noise of equation (II. lb has the form 

<Tjk(t)e- ihx (2.6) 

fceA<°°> 

where (/3fc) fcgA( oo) are independent <i-dimensional Brownian motions and ((Tfc) fceA (oo 
are d x d complex- valued matrices such that 

k ■ o k = (2.7) 

and 

kfc| 2 <oo. (2.8) 

Moreover, in order to obtain real-valued noise, we assume that 

aj = 0- k (2.9) 
for every k S A (oo) . Additionally, the vector-valued random field 

W(t,x)= VkMt)e~ ik - x 
fceA (oo) 

is, for every t > 0, partially isotropic if and only 

Wk\ = kflfcl (2-io) 

for all k G A (oo) and for every coordinate rotation R. 

Finally, in order to have measures with /i(L>(A)) = 1 we assume that 

J2 |£:|Vfc| 2 <°°, (2.11) 

since the values \k\ 2 con - espond to the eigenvalues of A. To summarise, we shall 
always assume that the noise <I2.6I > satisfies assumptions (I2.7t - (l2.11t . 
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2.2.2 The two-dimensional case 

The following result is well known. 

Lemma 2.8 Let jjl be an invariant measure of il.lt (d = 2) such that 

u j \\Du(x)\\ 2 dx < oo. 



r 



Then u G Vq and 

v ■ H \\Du(x)\\ 2 dx = ^ ^2 \ a k\ 



fceA(°°) 

/ 2 — ^ 2 2 

v - fjL I ||J5curlu(a;)|| dx = — I&I \ a k\ • 

It ^ — J 

Proof. 

Given fi, consider the (product) filtered probability space (Q, A,(At)t>o, P) 
supporting both a family of independent c?-dimensional Brownian motions 
(k, a) G A^ 00 ), and a non anticipating random variable no G Aq with law u. The 
corresponding strong solution u(t, x) of (ll.lt is a stationary process and satisfies, 
due to Ito formula, the balance relations 

\e p [ \\u(t,x)fdx + uE p [ [ \\Du(s, x)\\ 2 dx 
2 Jt Jo Jt 

= / \\ u o( x )\\ 2dx + 2 \ a k\ 2 -t 



\e p j \\cur\ u{t,x)\\ 2 dx + vE p [ [ \\D curl u(s, x)\\ 2 dx 
Jt Jo Jt 

l -E p ( ||curliio(s)|| 2 dz + i l^lVfcl 2 -*- 



2 

" 2 



fceA<°°) 

The result easily follows from stationarity. ■ 

Corollary 2.9 There exists a positive constant 6q, independent of v, such that 

0(p) > Oo 

for every invariant measure u G V of Hl.lt . 

Proof. The property > 9q follows from the definition of Q(p) and the two 
identities of the previous lemma, since 

j \\D 2 u(x)\\ 2 dx <C J ||£>curlu(x)|| 2 (fx 
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for a universal constant C > 0. ■ 

In the next theorem, when we say that M. C V x M + is a family of invariant 
measures of (II. lb , we clearly understand that each element (//, z/) G has the 
property that \i is an invariant measure for the Markov semigroup associated to 
equation (II . lb with viscosity equal to v. 

Theorem 2.10 In dimension d = 2, a family of invariant measures A4 C V x R + 
of il.lt cannot have any scaling law (in the sense of \2.2\ ). 

Remark 2.11 Under our assumptions on the noise, invariant measures of Hl.lt 
that belong to V certainly exist. In principle there could exist invariant measures 
for Hl.lt not belonging to V, but this has recently been excluded under very weak 
conditions on the driving noise (see M1M and the references therein). 

Remark 2.12 Consider equation rt7.il ) without the nonlinear term (called Stokes 
equations): 



in dimension d = 2, 3. Let M C V x R + be a family of invariant measures for 
it. Then the same results of the previous theorem hold true. The proof is the same. 
Alternatively, one may work componentwise in the Fourier modes and prove easily 
the claims. 

2.2.3 The three-dimensional case 

The lack of knowledge about the well posedness of the 3D stochastic Navier-Stokes 
equations has, among its consequences, the absence of the Markov property, and 
therefore of the usual notion of invariant measure. One may introduce several 
variants. Here we adopt the following concept. 

Consider the usual Galerkin approximations, recalled in Appendix B. The equa- 
tion with generic index n in this scheme defines a Markov process, with the Feller 
property, and has invariant measures, by the classical Krylov-Bogoliubov method: 
if X%(t) is its solution stalling from x and z/"' x is the law of X„(t) on H, by Ito 
formula it is easy to get a bound of the form (see for instance |5]]) 



which implies ([2| have been the first ones to use this elegant fast method) the 
necessary tightness in T of the time averaged measures 



If we choose the initial condition x = 0, then p,rf G V (in particular it is space 
homogeneous and partially isotropic), so there exist invariant measures in V for 



du 
~dt 



+ Vp = uAu+ <rkPk(t)e 



—ik-x 



fceA <oo) 
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the Galerkin equation. Denote by S n the set of all such invariant measures (thus 
S n c V). 

The constant C in the estimate above is also independent of n; it follows that 
the invariant measures of the class S n just constructed fulfill the bound 



< C. 



In fact it is possible to show that every element of S n has this property, J4j (if 
we do not want to use this property, it is sufficient to restrict the definition of S n 
in the sequel). These facts imply that U n S n is relatively compact in the weak 
topology of probability measures on H. We denote by V^^u) (the superscript G 
will remind us that we use the particular procedure of Galerkin approximations) 
the set of limit points of U n S n , precisely denned as follows: a probability measure 
fi on H belongs to 'PnsM if there is a sequence k n — > oo and elements fik n £ S kn 
such that fik n converges to fi in the weak topology of probability measures on 
H. The elements of the set V^ s (v) are space homogeneous and partially isotropic 
(these relations are stable under weak convergence). Furthermore, they have the 
other regularity properties required to belong to V: finite second moment in V 
comes from the previous estimates, fi(D(A)) = 1 from a regularity result of 0, 
see also |U, summarized in the following lemma. Therefore V^ s (iy) C V. 

Lemma 2.13 Given v > 0, there is a constant C v > (depending onv) such that 

^n(|^l5f) < C 

for every n and every invariant measure fj, n G S n . 

Given u G V , let S u be the tensor with L 2 (T) components 

S u = ^{Du + Du T ) 

(called stress tensor). The scalar field 

(S u (x)cur\u(x), curln(x)) 

describes the stretching of the vorticity field. If we set £ = curln, then formally we 
have 

g + ( u .V)f = i/Ae + 5'^ + i kxa k p k e~ ik - x . 
A formal application of Ito formula yields the inequality 

v ■ ji j \\Dcur\u(x)\\ 2 dx < fi J (S u (x)cur\u(x), curln(x)) dx 



+ 5 E i*r 



i2i i2 
-I \°k\ 
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for fj, G 'P§ s {v) (* n f act formally the identity). Along with the general results of 
the previous sections we would get 



x 



(S u (x)cur\u(x),cur\u(x))dx 



(2.12) 



This would be the final result of this section, having an interesting physical inter- 
pretation. However we are not able to prove it in this form. We analyze the status 
of this inequality by presenting some related rigorous results. They are of two dif- 
ferent natures: Corollary 12. 1 5l reformulates it for the coarse graining scheme given 
by Galerkin approximations; Corollary 12. 191 expresses the most natural statement 
directly for p G T^iy) DUt ^ re( l u i res an additional unproved regularity assump- 
tion. 

Lemma 2.14 Given fi G V^ s (u), an ^ ^n fe ^ S kn such that fj>k n converges to fi in 
the weak topology of probability measures on H, then 



\A- 



< lim/ 



\A- 



u 



The same is true for pjq- \\Dcur\u(x)\\ 2 dx in place of p 



\A- 



H 



Proof. Let {v?m} mg N G Cf,(H) be a sequence that converges monotonically 
increasing to \A-\ 2 H for every x G D(A), it is easy to construct it by cut-off and 
finite dimensional approximations). Since p(D(A)) = 1, by Beppo-Levi theorem 
lA<fm\ — > lA\A ■ |^]. Given g > 0, let mo be such that fi[(fm ] > lA\A ■ \ 2 H ] — e. 
Since /j n J(/3 mo ] — ► p[ip mo ] as k — > oo, eventually in k we thus have /inji^mj > 
p[\A ■ \ 2 H ] — 2s, and therefore also fi nk [\A ■ \ 2 H ] > p[\A ■ \ 2 H ] — 2e. This proves the 
first part of the lemma; the second one is similar. ■ 

Corollary 2.15 Let M CP x R+, with M u C V§ s (u), be a family with the K41 
scaling law, in the sense of deftnition U .3\ Then there exist Vq > and C > such 
that 



lim mfp nk 

k— >oo 



(S u (x)cur\u(x) , cur\u(x))dx 



T 



> c^ 2 (^ 



v • v 



-3/2 



for every v G (0, vq), every p G M. v and every sequence p n , k G S kn such that pk n 
converges to p, in the weak topology of probability measures on H. 

Proof. From the previous section we know that 

p J \\D 2 u(x)\\ 2 dx > e 3/2 (p,v) ■ v~ 5/2 . 

Since 



(Af,g) H = (cur\f,cur\g) H 



(2.13) 
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for every /, g € D(A), we have 

H j \\Dcur\u(x)\\ 2 dx > Ce 3/2 (fi,v) ■ v~ 5 / 2 
for a suitable universal constant C > 0. From the previous lemma we have 

liminfy nfc / \\Dcur\u(x)\\ 2 dx > Ce 3/2 (fi, v) ■ v^ /2 . 
Thus the claim of the corollary will follow from the inequality 

v ■ fj, nk J \\Dcur\u(x)\\ 2 dx < [i nk J (S u (x)cur\u(x), curlu(x)) dx (2.14) 

Let us sketch the proof of this inequality (see |4) for more details). Consider the 
Galerkin approximations 

du {n) +\uAu^ +^B[u^ n \u^y^dt= £ °kd(3 k e- ik - x 

fceA<») 

described in Appendix B. From Ito formula for (Au( n \t), u^ n \t)) H we get 
(Au {n \t),u {n \t)^ H + / 2^ M W,i/iii<"' + ir^B(u {n \u^^ H ds 

= (A«W(0),uW(0)) ff + MT + i £ |fc|V fc | 2 

fceAg°> 

where M t n is a square integrable martingale. We have 

since 7r( n ) is selfadjoint and commutes with A. Besides d2.13t we also have 
(Af,B(g,g)) H = (curl/, (g • V)curl 5 + S g cur\g) H 

hence 

(Af,B(fJ)) H = (cur\f,S f cur\f) H 
for every f,g£ D(A). Therefore we have 



cu 



< 



r\u {n) (t) 2 H + J Uu | DcurluW | ^ + (curlu^ , 5 uW curlu (n) ^ ds 



This implies (12.141) and the proof is complete. 
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Remark 2.16 We cannot conclude \2.12\ from the previous corollary without fur- 
ther (unproved) assumptions on p or {p nk }- This could be just a technical point 
due to the present lack of better regularity estimates for the 3D Navier-Stokes equa- 
tions, or it could be a facet of a deeper phenomenon. Let us explain it with a 
cartoon argument. First recall that it is easy to construct, say on the torus T, a 
sequence {f n } of functions converging a.s. to zero, but with f n dx = 1 (or even 
Jj- fn dx — > oo ): just take the modifiers of a Dirac delta distribution; if we like, the 
example can be modified so that f n tend to develop singularities on a dense zero 
measure set in T, but the a.s. limit is still zero. Thus we see that for the limit mea- 
sure \i we could have a small value of fi [f T (S u (x)cur\u(x), curlu(x)) dx] even 
if some coarse graining procedure, here represented by the Galerkin approxima- 
tions, could give us a large value of p Uk [f T (S u (x)cur\u(x), curlu(x)) dx]. Such 
arguments arise the question of the physical meaning of the true Navier-Stokes 
equations and possibly of its coarse graining approximations; this is not our aim, 
but we wanted to say that the previous corollary may be considered perhaps as a 
result of possible physical interest in itself, even if we cannot rewrite it in the form 

Lemma 2.17 Given fi £ ^ > Ns( zy )' an d ever y sequence p Uk G S kn such that /j,^ 
converges to fi in the weak topology of probability measures on H, we also have 
p nk — > n weakly on [W 1,3 (T)] . 

Proof. From the lemma above, {^n k } is bounded in probability on D(A): 

^ nk (\Ax\ H >R)= » nk (\Ax\f > < R-V^ n (\A-\ 2 f 

C 

< 



R 2/3- 

The embedding of D(A) into [VF 1,3 (T)] 3 is compact: recall that Sobolev embed- 
ding theorem gives us 

W 2,2 c W P,2f3=T for every /3 G (1,2), and the embedding of 
in W '2/3-1 is compact; choose then (3 = 3/2. Therefore {f^n k } is tight in 
[W^iT)] . Easily we deduce that it converges weakly to fi also in [W /1,3 (T)] . 
■ 

Corollary 2.18 If fi G Vg s (v) is 

the weak limit (in H and thus in [M /1 ' 3 (T)] ) of 

a sequence /x„ fe G S kn such that 



|2+£ 

ly 



< C 



for some e, C > 0, then 

|2 



v-p\ ||£>u(»|| 2 (i:r = ^ ^ \a k 

r fceA<°°> 
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If in addition 



then 



Mn fe 



|3+e 
W 



< c 



i J ||jDcurlit(z)|| dx < fx J (S u (x)cur\u(x), curln(x)) dx 



+ 



- 2 £ i*i 

fcgA (oo) 



2| i2 
0~k\ ■ 



Proof. It is sufficient to apply repeatedly the following fact: if /i n — > p weakly 



< C, then fi n {p] — > fi[tp]. This 



in a Polish space X, ip G C(X) and /U n M 1+e 
fact is well know but we provide the proof for completeness. Let Y n and Y be 
r.v.'s with law [i n and resp., with values in X, such that Y n — > V a.s. in X. 
Then /i n [</?] = E[tp(Y n )], p[tp] = E[tp(Y)], so by Vitali convergence theorem it is 
sufficient to prove that <p(Y n ) is uniformly integrable. We have 

E[<p(Y n )l v{ y n) > x } < {EMY n f]flvp{v{Y n ) > A) 1 / 9 < CX- 5 . 
Thus the uniform integrability is proved and the proof is complete. ■ 

Corollary 2.19 Let M CP x R+, with M u C V§ s (v), be a family with the K41 
scaling law, in the sense of definition \1.3\ Assume that every \x in M is the weak 
limit of a sequence fi n , G S kn such that 



fi>n h 



|3+e 



< c 



for some e, C > 0. Then there exists v$ > and C > such that \2.12l holds for 
every v G (0, Vq) and every \i G M. v . 

Remark 2.20 If K41 scaling law holds then vortex stretching must be intense. 
Heuristically, no geometrical depletion of such stretching may occur ( in contrast to 
the 2D case where the stretching term is zero because curlu(x) is aligned with the 
eigenvector of eigenvalue zero of S u (x)): indeed, if we extrapolate the behaviour 



E 



\Du\ 



^ - as Du ~ — t=, curlu ~ -4=, then we get E\S u cur\u ■ curlul ~ 
if there is no help from the geometry. Another way to explain this idea is the fol- 
lowing sort of generalised Holder inequality. 

Corollary 2.21 Let M CP x R +) with M v C V§ s (u), be a family with the K41 
scaling law, fulfilling the assumptions of corollary \2.19\ Then there exists uq > 
and C > such that 



\ 1/2 

H I \\Du\\ 2 dx ) <C[p 
for every v G (0, vq) and every p, G M. v . 



\S u zw\u • curlu| 



'dx 



T 



1/3 
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Proof. From the previous corollary and the definition of e(/i, v) we have 



(i 



\S u cw\u ■ curl-u|| 2 (ix 



r 



V3 > (Ce 3 / 2 ^,^)-^- 3 / 2 ) 173 



1/2 

\\Du\\ dx ) 

T 



The proof is complete. ■ 

3 Necessary and sufficient conditions for K41 

We continue with the notations and concepts just introduced in the last section on 
the 3D case. 

The result of this section can be formulated for definition 1 1.31 but the presence 
of the factor e(/x, v)- 1 / 4 in the definition of rj(p,, v) makes some statements much 
less direct. So, having in mind the exploratory character of these equivalent con- 
ditions, we prefer to adopt a simplified form of our definition of the K41 scaling 
law. 

Definition 3.1 We say that a scaling law of K41 type holds true for a set M C 
V x M + if there exist vq > 0, C > c > 0, Co > 0, and a monotone function 
Rq: (0, vq] — > with Rq(v) > Cq and lim^^o Roiy) = +oo, such that the 
bound 

c-r 2 / 3 < S%{r) <C -r 2 / 3 (3.1) 
holds for every pair (fi, v) G M. and every r such that v G (0, vq\ and 

C u 3 / 4 < r < ^ 4 Rq(v). 

Recalling that rj(fj,, u) = v 3 / 4 e(n, v)~ l l 4 , we see that this definition is equiva- 
lent to ll.3l if there exist e\ > eo > such that 

eo < < ei 

for all (fj,, v) G M.. Unfortunately, in 3D only the upper bound can be proven. 
However, this could be just a technical problem due to the fact that we can only 
use weak solutions (for slightly more regular solutions Corollary 12. 18l implies that 
e(fj,, v) would be bounded from above and below). 

Consider the auxiliary stochastic Navier-Stokes equations 

-£(t, x) + (u(t, x) ■ V)u(t, x) + Vp(t, x) = 9Au(t, x) + V a k $ k (t)e- ik - x 

(3.2) 
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on the torus [0, L] 3 with div u = and periodic boundary conditions (the set A^ o) 
is defined in (IB.lt ). As we shall see below (see the next section and lemma lB~lT i. 
we obtain this equation when we perform the following scaling transformation on 
the solutions u of the original equation WAY . 

u(t,x) = L l ^u{L~ 2 ^t, L~ l x) 

(and a suitably defined p(t, x)). The value of v under this transformation is 

v = vL^I\ 

This scaling transformation has been introduced in the mathematical-physics liter- 
ature, see Q. What makes it special is that no coefficient depending on the scale 
parameter appears in front of the noise, so the energy input per unit of time and 
space is the same for every L. Heuristically, if we believe in a cascade picture of 
the energy (without essential inverse cascade), this invariance of the energy input 
should imply that the small scale properties of (11.1b and (13.21) are the same, namely 
that they are invariant under this transformation; this should lead to the K41 scaling 
law. 

Similarly to the case L = 1, we may introduce the (non empty) set V^ S (D, L) 
of limit points of the (homogeneous and isotropic) invariant measures of the corre- 
sponding Galerkin approximations. 

Let us denote by "P^s the set of all pairs (p, v) such that p G V§ s (u). Similarly, 
let us denote by V§ s the set of all triples (p, u, L) such that p 6 V^ s (i>, L). 

3.1 Basic equivalent condition 

The following condition seems interesting since it looks rather qualitative, in con- 
trast to Definition 13.11 and shows that the exponent 2/3 arises from the scaling 
properties of the stochastic Navier-Stokes equations. 

Let us introduce the notation Vl for the set of probability measures analogous 
to V, but on the torus [0, L] 3 . Denote by V. x the set of all triples (p., v, L) 
such that L) £ and fi G Vl- In the next definition and later on we use the 

u(e) — u(0)\\ 2 when fi G Vl (and other similar mean values): this 



notation /j, 
means 



|u(e)-u(0)|r = / ||u(e)-u(0)H , 



1 1 

where Hl has been introduced in section [LTl 

Definition 3.2 We call admissible region a set D C of the following form: 

D = \(p,L) G K+; v G (0,u o ),L> R (9)} 

where Dq > and Rq: (0, vq] — > [l,oo) is a strictly decreasing function with 
Ro(v) — > cxo as v —> 0. 
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An admissible region is depicted in the left-hand side of Figure Qbelow. 

Condition A A subset M. C V. x is said to satisfy Condition\A\if there exist 
an admissible region D C and two constants C > c > such that 



c < \i 



|«(e)-«(0)|| : 



< C (3.3) 



for every (/i, D,L) £ M with (P, L) G Z). 



Proposition 3.3 The set Vg satisfies Condition^if and only if the set VS S has a 
scaling law ofK41 type, in the sense ofDeftnition \3.1\ 

Proof. Given R > 0, consider the mapping Sr : TLr — ► H defined by 

(S R u)(x) = R l/3 u(Rx) . (3.4) 

This mapping induces a mapping S from V x to? x R + by 

SQi,u,r) = (5i/i,z>f- 4 / 3 ) . (3.5) 

It follows immediately from Theorem IB . 21 that one has 

Vg=S(Vg), and pg = S~\vg). (3.6) 

Furthermore, it follows immediately from the above definitions that if (//, v) = 
S(jl, u, f ), then 

S$(r) = r 2/3 f \\u(e) - u(0)\\ 2 dft(u) . (3.7) 

It therefore follows that, in order to prove the equivalence between Condition |X] 
and K41, it suffices to show that the domains of validity of eq. l3.3l and of eg. 13. II 
are the same (with possibly different constants and functions Rq and Rq), provided 
that {y, r) and (p, f) are related by 

v = ur~ 4/3 , r = r" 1 . (3.8) 

We denote by K : (y, r) \— > (P, f) the above map. 

Condition lAlimplies K41. The domain of validity of eq. r3.3l is given by 

v < z>o , r> Rq(u) . (3.9) 
Under the map K~ l , this becomes 

r > ( ^) 3/4 = Co v 3/A , - > R (vr~^) . (3.10) 
\vq/ r 

Both domains are shown in Fig.^ 
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Parameter domain for Condition A Image of the previous domain 



Figure 1 : Effect of K 1 on an admissible domain 

Denning the strictly decreasing function F(x) = x~ 3 ^Ro(x), the second con- 
dition of eg . 13 .101 is of course equivalent to 

„-3/4 > F(vr~ 4/3 ). (3.11) 

This condition (as can be inferred from the Fig. can only be satisfied simulta- 
neously with the first condition in eq. 13 .101 if v < uq = F(i>o) -4 ^ 3 - On (0, vq] this 
domain, eq. 13.1 H is equivalent to 

r<( jrT ^y / ^^R Q (,), (3.12) 

where Ro(x) = (F~ 1 (x~ 3 / 4 )) -3 ^ 4 . Additionally i?o is well-defined on (0, uq\ and 
that it is greater than Co on this domain. Furthermore, since F is decreasing, Rq 
is strictly decreasing and it is easy to check that lim^o Ro(%) = oo because the 
same property holds for F. 

K41 implies Condition EJ The domain of validity of K41 is given by 

v<vq, rv~ z l A e [C , Ro(is)] . (3.13) 
Under the map K, this becomes 

£r" 4/3 < , z>" 3/4 G [C , i? (^- 4 / 3 )] . (3.14) 
The second condition can be rewritten as 

D£ [C(z>f- 4 / 3 ),£ ] , (3.15) 
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Figure 2: Effect of K on a domain of the type < I3 . 1 3 1 

where we defined uq = C 4 ^ 3 and G(x) = Rq(x)~ a ^ 3 . Both of these domains are 
shown in Figure |2] 

We can rewrite as above the condition v > G(Pf~ 4 / 3 ) as 

r *(o=m) sR » (f) - (3 - 16) 

Again, it is an easy exercise to show that Rq as denned above is monotone and 
satisfies lirn^o Ro( x ) — °°- The only points that remain to be clarified are: 

a. We haven't taken the first equation in eq. l3.14l into account. 

b. The domain of definition of Rq may not extend to Po- 

Both problems can be solved at once by simply choosing a smaller value for i> ■ 

3.2 Necessary and sufficient conditions in terms of high and low modes 

Although Condition |A] contains only bounds (at finite distance points) and not 
scaling exponents (with small distant points), and thus in principle it represents 
a progress in the direction of analysis of K41, it still looks difficult to verify or 
disprove it for Navier-Stokes equations, since it is rather unusual to work with the 
difference of a solution at two points. This is the main motivation for the following 
new necessary and sufficient conditions. 

Looking at them on the other direction, as necessary conditions for K41, they 
declare that under K41 the energy of high modes is bounded and the enstrophy of 
low modes is bounded, an information with a certain physical content. 
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In this section, for notational simplicity, we drop the tildes in our notation. 
Recall that an admissible region is defined by 

D = {(u,L)€R%;v€(p,vo),L>Ro(u)} , 
and that Condition lAlrequires 



c < \i 



|u(e)-u(0)||' 



< C 



for every (/i, v, L) with (y, L) G D. 

We start with a preparatory lemma which depends on the scaling properties of 
stochastic Navier-Stokes equation in an essential way. This is the only point in this 
section where specific informations about the measures are used. 

Lemma 3.4 If V^ s satisfies Condition A then there exist constants C > d > 
and an admissible region D' such that 



|u(Ae) - u(0)||' 



dX < C' 



for every (fi, u, L) G V^ s with (u, L) G D'. The sum ^ e is extended to all coordi- 
nate unitary vectors. We simply have C = (l.5 2//3 d) • C, d = (0.5 2 / 3 <i) • c, D' 
defined by 0.5 4 / 3 • Vq and 1.5Rq (l.5 _4//3 ^), where uq and Rq(v) define D. 

Proof. Given A G [|,|] and (fi,u,L) G fi§ s , namely n G V^ s {u,L), 
consider the measure n\ that corresponds to n under the transformation u \— > 
A _1 / 3 u(A.) used in the previous section, having the property 



|u(Ae) -«C0)|| ! 



A 2/3 



/'A 



|u(e)-u(0)|| : 



By TheoremElwe know that fix G 7^s(^A~ 4 / 3 , L/X), hence (/x A , i^A" 4 / 3 , L/X) G 
"P^s- Thus Condition lAlimplies 



C < ^A 



it e 



u(o)ir 



< c 



if z^A -4 / 3 < i/o and L/A > i?o(^A -4 / 3 ). The first condition is true if v < 
0.5 4 / 3 z/o- The second one if L > 1.5i2o(l.5 -4 / 3 z/). The proof can now be easily 
completed. ■ 

Let us use some Fourier analysis on the torus Tl = [0, L] d (see also Appendix 
iBl . Every u G Hl is given by 



u{x) 



ik-x 



u{k) 
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where 



u{k) := L~ A 




e u(x)dx 



and we have Parseval identity 




u{x)\\ 2 dx= E R*0f 
keA^ 



We introduce another condition which requires the sum of the enstrophy of low 
modes and energy of high modes to be finite and bounded away from zero. 

Condition B A subset J\A C V. x is said to satisfy Condition\E\if there exist 
an admissible region D C M. 2 , and two constants C > c > such that 



for every (fx, z/, L) £ V^ s such that (u, L) £ D. 

With this definition, we may establish a first basic theorem as a corollary of the 
previous lemma. 

Theorem 3.5 Condition \A\inwlies Condition \B\ 

Remark 3.6 We understand that constants and admissible regions involved in Con- 
ditions$\and\E\are not necessarily the same. 

Proof. For every u £ Hl we have 



c < 



E wkfvlwmw 2 } + E v[\\m\\ 2 } ^ c 





E k 



,ik-Xe 



1 Rfc)|| 2 



and thus, for every fi £ V^ s (^, L) we have 




e 
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But there exist universal constants C > d > such that 

2 







I: 









dX < C'(\\k\\ 2 A 1). 



Therefore, the quantities 



and 



£ / /i[||u(Ae)-«(0)|| : 

e ''I 

E (iifcfA^fiiswir 



fceA^ 

are "equivalent", up to universal constants. This proves the claim. ■ 

We have at least a partial converse of the previous result if we require that in 
the admissible region the enstrophy of high modes is by itself bounded away from 
zero. Then we introduce the following condition: 

Condition C A subset M C V. x IR^_ is said to satisfy Condition\C\if there exist 
an admissible region D C and two constants C > c > such that 



c< E \\k\ 



fceA 



(oo) 



< E iniv 

fceA 1 



Ki 



m\\ 2 }+ E 

fc£A (oo) 



< c 



fcll>l 



||fc||<l/2 

for every (fi, u, L) £ "P^ s smc/j (z/, L) G Z). 

Note that Condition |C] implies directly Condition |B] What is more interesting 
is the following: 

Proposition 3.7 Condition \C\implies Condition^ 
Proof. We have 



E 



1 



< C7f HArlr A 1 



for every k. Moreover if \\k\\ < 1/2 we have 



< 



E 



e ik-e _ 1 



for some constant c > 0. The claim then follows from the next lemma and the 
following inequality 



E Ek"- 1 



fceA 



(oo) 



I/' 



* E (E 

A,-eA ( r 



e ik-e _ 1 



/' 



\u(kW 



||fc||<l/2 



E 

fcGA 



(oo) 
L 

fc||<l/2 
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Lemma 3.8 V^ s satisfies Condition A if and only if it satisfies the following Con 
dition A': there exist C > c > 0, and an admissible region D such that 



for every (fi, u, L) G V^ s such that (y, L) G D. 

Proof. From previous computations, we know that for every fi G V^ s (u,L) 
we have 



This proves the claim. ■ 
4 A random eddy model 

We now exhibit a model having the property stated in the conjecture, and other 
heuristically meaningful properties for a turbulent velocity field. The model is 
mathematically rigorous but it is not derived from the Navier-Stokes equations, it 
is just a cartoon of what we believe to resemble the turbulent 3D field given by 
the Navier-Stokes equations. Therefore the only merit of the following result is to 
show that there exists a field with the property stated in the conjecture, and such 
a field is not just an artificial example but it is strongly inspired by numerical and 
physical observations of turbulent fluids. 

For simplicity we work in the full three-dimensional space M 3 , instead of the 
torus T. 

The model should be thought of as a random collection of vortex filaments, 
i.e. concentrations of vorticity around one-dimensional continuous curves. The 
filaments will be of various kind, from very elongated ones, whose existence is 
well documented in numerical observations of fully developed turbulence, to other 
more "eddy-like" and symmetric. 

The basic ingredient of the construction is a vortex filament of length T, thick- 
ness I and core velocity U, which is stochastically modelled around a "Brownian" 
core: consider a 3d-Brownian motion {X t }te[o,T] starting from a point Xq. This is 
the backbone of the vortex filament whose vorticity field is given by 





2 




(4.1) 
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where odX denote Stratonovich integration. The letter t, that sometimes we shall 
also call time, is not physical time but just the parameter of the curve. We as- 
sume that gt(x) = g(x/£) for a radially symmetric measurable bounded (smooth) 
function g with compact support in the ball B(0, 1) (the unit ball in 3d Euclidean 
space). Heuristically ^ s i ng ie(x) is an average of the "directions" dX t for points X t 
in the ball B(x, £). The various parameters U, £, T have to be thought of as giving 
the "typical" magnitudes of the respective properties. It should be noted that £ is 
not a "real" vorticity field (since in this model its divergence is not zero) but should 
be understood as providing the contribution to the fuild vorticity coming from the 
eddies. 

The velocity field u is generated from £ according to the Biot-Savart relation 

U f T 

Usingle(x) =p J K-1.{X - X t ) A 0dX t (4.2) 

where the vector kernel is defined as 

1 f x — y 

fci(x) = — Qi(yh Tndy. (4.3) 

47r JB(0,£) k — 2/1 

We want to describe a random superposition of infinitely many independent Brow- 
nian vortex filaments, uniformly distributed in space, each of which will be asso- 
ciated with intensity-thickness-length parameters (U, £, T) "randomly drawn" ac- 
cording to a measure 7. The total vorticity of the fluid is the sum of the vorticities 
of the single filaments, so, by linearity of the relation vorticity-velocity, the total 
velocity field will be the sum of the velocity fields of the single filaments. 

The correct mathematical implementation of this heuristic picture is given by 
the construction of a Poisson random measure on a suitable space. 

Let S be the metric space 

5 = {(U,£,T,X) G M 3 x C([0,1];M 3 ) : < £ < VT < 1} 

with its Borel <7-field £>(E). Let (Q, A, P) be a probability space, with expectation 
denoted by E, and let [i^, u £ fl, be a Poisson random measure on £>(H), with 
intensity v (a cr-finite measure on /3(H)) given by 

du(U, £, T, X) = dj(U, £, T)dW{X). 

for 7 a cr-finite measure on the Borel sets of {(U, £, T) G : < £ < \Pf < 1} 
and dW(X) the cr-finite measure defined by 



f i>(X)aw(X) 

JC([0,1];R 3 ) 



TP(X)dW X0 (X) 
C([0,1];R 3 ) 



dxr 



for any integrable test function ip : C([0, -> R. Here dW Xo (X) is the 

Wiener measure on C([0, 1],M 3 ) stalling at xq and dxo is the Lebesgue measure 
on M 3 . Heuristically the measure W describes a Brownian path starting from an 
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uniformly distributed point in all space. The assumptions on 7 will be specified at 
due time. 

The random measure /x w is uniquely determined by its characteristic function 



Eexp / ipiQuuidQ = exp / (e 



MO 



for any bounded measurable function 92 on E with support in a set of finite u- 
measure. In particular, for example, the first two moments of \i read 



E / tfKOKdO 



<p(Qu(dO 



and 



E 



tp(Qv(dQ 



+ 



<p 2 (Qv(dO. 



Given the Poisson random measure jjl we can introduce our random velocity 
field as 

u(x) = u!? ingle (x)MdO = M(«kngto(a?)) • (4-4) 
for any x G M 3 , where C = (U, £, T, X) and 

C " Single (?) ■= JI [ ~ X t) A dX t 

can be shown to be a well defined //-measurable function. 

In plain words, given uj G 0, the point measure fj,^ specifies the parameters 
and locations of infinitely many filaments: formally 



(4.5) 



for a sequence of i.i.d. random points {C"} distributed in S according to u (this 
fact is not rigorous since v is only cr-finite, but can be justified by a localisation 
procedure). Since the total velocity at a given point x G M 3 should be the sum of 
the contributions from each single filament, i.e. in heuristic terms 



u{x) = XXiie( x ) 



(4.6) 



this justifies, physically, the above formula. 

To end the construction of the model it remains to choose a suitable measure 
7 for the distribution of the parameters. Lacking physically motivated choices of 
7 we resorted in Q to show that it is possible to fix 7 in such a way to recover 
statistics which corresponds to multifractal scaling of the velocity increments, for 
any possible choice of the multifractal spectrum. In this way we showed how 
to build a random field with prescribed multifractal spectrum which also possess 
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some geometric properties of real turbulent fields. In particular we can choose 7 to 
recover K41 behaviour of the velocity increments. In the following we will fix this 
particular choice and show that, for our random field 

E\Du(0)\ 2 ~ ??" 4/3 E|D 2 it(0)| 2 ~ 77 10/3 (47) 

where rj, in the context of this section will be a UV cutoff scale for the vortex 
model, i.e. we will not allow vortices with thickness t smaller than 77 which physi- 
cally models the "viscous" (or Kolmogorov) scale which determines the lower end 
of the inertial range. The scaling (14.7b implies that 9 ~ 77 for 77 — > 0. 
So we stipulate that 

d 7 (U, £, T) = 5 el/:i (U)5AT)rH mvA) de (4.8) 

where we ignore a possible constant prefactor which will not play any role in our 
discussion. This choice of 7 corresponds to force the vortex filaments with thick- 
ness I to have length proportional to £ 2 and to have typical velocity of the order 
of i 1 / 3 , the "density" £ -4 is chosen to roughly have "space-filling" vortices at all 
scales. Moreover vortices can have thickness going from the small scale 77 to a 
large "integral" scale of order 1. 

Let us state the result. For technical reasons we will assume that there exists 
positive constants c, C, A and u G B(0, 1) such that the following bounds on K, 
holds 

cl x eB(u,\) < \DK x (x)\ < Cl xeB (.o,D (4-9) 

and 

cl x&B (u,X) < \D 2 Jd(x)\ < Cl xeB( p tl) . (4.10) 

for any x G B(0, 1). A sufficient condition for the upper-bounds is that g is 
bounded. 

Proposition 4.1 With the above definitions, we have 

E\Du(0)\ 2 x r/" 4/3 K\D 2 u(0)\ 2 X r/~ 10/3 (4.11) 

as 7] — > 0. 

Proof. By a small abuse of notation, we have 

t2 r T 



E\Du(0)\ 2 = v[\Du siagle (0)\ 2 ] =u 



\D)C e (0-X t )\ 2 dt 



since K-Du S mgi e (0)) = being n s i ng i e a v-lto integral and where we used the 
energy-identity for the Ito integral. Note that \DKLi(x)\ = \DJC\(x/£)\ and that, 
when x B(0, 1) we have 

\Dld(x)\ < C\x\~ 2 , |£> 2 /Ci(x)| < C7|x|~ 3 (4.12) 
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by direct estimation from the Biot-Savart formula. 

Now we use the Lemma l4~2l below together with the bounds < I4.9I > and (14.121) to 

get 



E|£(0)| 2 x 7 [[/ 2 r 1 T] =7 



^5/3 



where we used the fact that, under 7, U = I 1 ' 3 and T = I 2 . Then easily we 
conclude that 

e|zmo)| 2 - ?r 4/3 - 



Analogously, from the bounds (14. 101) and (14.121) we have that 

cr l i x&B(uM) < \D 2 ic e {x)\ < cr 1 [i x€B (o,£) + i^B(o,d^r 3 ] 

< C'l' 1 [l x eB(0,£) + lx&B(0,£)\x /£\~ 2 \ 

so using again Lemma |4~21 we can obtain that 



E\D 2 u(0)\ 2 x -y[U 2 r 3 T] 



7 



£-1/3 



-10/3 



ending the proof. 



Lemma 4.2 We have the estimate 

w' 



1 







x t eB(o, 



and, if 
we have 



i(x)\ < C{l xeB (o,e) + lx#B(.ojb\x/ 



\ip t (x - X t )\ 2 dt 



< tT 



(4.13) 



(4.14) 



(4.15) 



Proof. These results are particular cases of more general bounds proved in Q : 
the first is proved in Lemma 14 of the reference, while eq. d4.15t is proved in 
Lemma 3: the proof refers to the particular case in which (pi = Kg but it is easy to 
see that a sufficient condition is given by eq. (!4.14l) . ■ 



Appendix A Mollification of measures 

Some computations of the paper with Taylor formula require more regularity than 
that of typical fields under \i € V. For this reason we introduce mollifications of 
measures \i G V. Let us remark that this technical effort is useless if the noise is 
more regular, since one can prove more regularity of the typical elements under 
H G V. 

Let ip : E — > R be a smooth function with compact support, symmetric, 
non negative, strictly positive at zero, with f Rd cp(\\x\\)dx = 1. Set p e (x) = 
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e _d V?(||a;/e||), so J Rd ip £ (x)dx = 1; {ip £ } £>0 is a family of usual smooth mol- 
lifiers. For every u G H set 

u e (x) = / ip £ (x - y)u{y)dy. 

Given [i G Vq, the mapping u ^ u £ in H induces an image measure fi £ G Vq 
which is in fact supported on smooth fields. 

Lemma A.l If fi G V, then fi £ G V. 

Proof. We have 

u £ (x-a)= cp £ (x - a - y)u(y)dy v ~= +a ip £ (x - y')u{y' - a)dy' 



£ 



Ve{x - y')u{y')dy' 



where the last equality is understood in law under fi, and it holds true as processes 
in x. Hence u £ {- — a) = u £ (-). This means 

f(u £ {- - a))d^{u) = I f{u £ )d\i{u) 

H JH 

for bounded continuous /'s, and therefore 

J f(u(- - a))dfj, s (u) = [ f(u)dfi £ (u) 



11 



so the space homogeneity of fi £ is proved. 
Similarly, we have 

u £ (Rx) = (p £ (R(x - R~ 1 y))u{y)dy = ip £ (x - R~ l y)u(y)dy 



form the symmetry of tp £ , hence 



u £ (Rx) y = v I tp e (x - y')u{Ry')dy' 



c 



Ve{x - y')Ru(y')dy' 



hence 



and finally 



f(u £ (R-))dfi(u) = / f{Ru E {-))dn(u) 

H JH 



The proof is complete. 



(R-))dfi £ (u) = [ f{Ru{-))d^ £ {u). 

JH 
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Lemma A.2 For every \i EV, if 




\\Du(x)\\ dxdfi(u) < oo, / / \\D u(x)\\ dxd/j,(u) < oo, 



then 



and 



/ \\u(re) — u(0)\\ 2 dfj,(u) < oo 
Jh 



e^O 



lim / \\Du(0) \\ 2 dn £ (u) 



H 



\\Du(x)\\ 2 dx dfi(u) 



H JT 



lim / \\D 2 u(0)\\ 2 dfi £ (u) = [ [ \\D 2 u{x)\\ 2 dx dfi(u) 
^°Jh JhJt 

lim / \\u(re) -u(0)\\ 2 dn £ (u) = / ||u(re) - u(0)\\ 2 dfi(u) . 
£ ^°Jh Jh 



Proof. There exists C > such that 

2 



\\Du £ (x)\\ z dx < C / \\Du{x)\\ A dx 

T Jt 

for every u G D{A); and J T \\Du £ (x) \\ 2 dx — > f T \\Du(x) \\ 2 dx as e — > for every 
it G £>G4). Hence, by Lebesgue theorem, 



lim 



\\Du(x)\\ 2 dx 



T 



dn e (u) 



\\Du(x)\\ dxdfx(u). 



H JT 



But /i £ is space homogeneous, hence 




IT 



||Du(x)|| 2 da: 



dfi e (u) = I \\Du(0)\\'d^(u) 

H 



This proves the first claim. The proof of the second one is entirely similar. For the 
third one, we have 



\u £ (x + re) — u £ {x) 



r / Du £ (x + ae)eda 
Jo 

<r 2 \\Du £ (x + ae)\\ 2 do- 
Jo 



for every u G D(A), hence 



u £ (x + re) — u £ {x)\\ dx < r / \\Du £ (x + ere) || dxdo 

Jo Jt 

= r 2 J J \\Du £ (x)\\ 2 dxdo 
<Cr 2 J \\Du(x)\\ 2 dx. 
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Therefore, again by Lebesgue theorem, 

\u(x + re) — u(x)\\ 2 dx dfi £ (u) 
\u(x + re) — u(x)\\ 2 dx dfi(u) . 



lim 



H JT 



H JT 

The third claim follows now from the space homogeneity of both \i e and \i. ■ 

We are now in the position to prove a quantitative consequence of isotropy, that 
we shall use in the sequel. In the next statement we understand that both terms in 
the equality are either finite and equal, or both infinite. 

Lemma A.3 For every [i £ V and every coordinate unitary vector e we have 



H JT 



\Du{x)\\ 2 dx d/j,(u) = d 



H JT 



\\Du(x) ■ e|| dx dyi{u). 



For n e , we have the same identity and also 



\Du 



H 



dfj, e (u) = d I \\Du{0) ■ e|| d^ £ {u). 



H 



Proof. Step 1. Denote by ei,...,ed the coordinate unitary vectors. For u G 
D(A) we have 



\Du(x) 



E 



dxj 



(x) 



and thus 



\Du(x) 



\Du(x) • e j- 1 1 2 = 



\\Du(x) ■ ej || 2 . 



duj 
dxj 



x) 



Therefore 



and 



\Du 



ii 



]T I \\Du{Q)-ejfd^u) 



j j ||£>ii(a:)|| 2 da;d/z(u) = 

J H JT 



H JT 



\Du{x) ■ ej\\ 2 dx djjb{u). 



It is then sufficient to prove that all terms of the sums on the right-hand-sides are 
equal, in order to prove the first and last claim of the lemma; we shall prove this 
below in steps 2 and 3. Finally, the first assertion for /i e is a particular case of the 
first claim of the lemma (jj, s is an element of V). 
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Step 2. Now, given j = 1, d, take a rotation R as in the definition of V such 
that Rei = &j. Given N > 0, 



Du(0) • ej || 2 A Njdfi £ (u) 

of ( r ~ 2 \\ u ( re i) ~ u (°)\\ 2 A N ) d ^( u "> 
of ( r ~ 2 \\ u ( Rre i) ~ U ( R0 )\\ 2 A N ) d ^( u ) 
lim / (V _2 ||u(rei) -u(0)|| 2 A JvW e (u) 
Du(0) -eif AN\d[i E (u). 



u 

lim 



lim 

r 



By monotone convergence in N, we get that j H \\Du(0) ■ ej\\ 2 dfi e (u) is indepen- 
dent of j. This proves one of the claims. 

Step 3. From the previous step and homogeneity we have that j H j r \\Du(x) • ej \\ 2 dx dfi e (u) 
is also independent of j. Arguing as in the proof of the previous lemma, this inte- 
gral converges to j H j T \\Du(x) ■ ej\\ 2 dx dfi(u), which is therefore also indepen- 
dent of j. The proof is complete. ■ 

Appendix B Scaling theorems 

The torus, Tl = [0, L] d , the energy space Hl with norm , the spaces Vl, 
D{Al), T>l and the Stokes operator Al on Tl have been already introduced in 
section [TTT1 We define 

Af> = Le : \k\ 2 >0 J, (B.l) 

and, for the purpose of Galerkin approximations, we introduce also 

A^ = jfce^z d . <\k\ 2 < {^-n 2 

so that a£° } = U n A^ n) . In particular, A (oo) = A<f*>\ 
B.l Scaling theorem for Galerkin approximations 

Let V' L be the dual of Vl', with proper identifications we have Vl C Hl C V' l 
with continuous injections. Let Bl(-, ■) '■ Vl x Vl — > V' L be the bilinear operator 
defined for all u,v,w £ T>l as 

x^ 1 f dv, , \ - 



(w,B l (u,v)) Hl = 2^ j2 J u iQ^ w J dx = 2^ {I -u(h))v(l)-w(k). (B.2) 

i,j=l ^ L % h+l=k 
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Given L > 0, v > and 9 > 0, consider (formally) the equation in Hl 

du+ [vA L u + B L (u,u)]dt = 6 ^ of e~ ik ' x , 



where (3% = (3 L k and o f = cr Lfe , and (/3 fc ) fcgA(0 o) and (d fc ) fcgA ( o) have been intro- 
duced in Section l2.2.1l and are subject to the assumptions imposed therein, so that 
the random fields 

W L n \t,x) = £ o L k (l L k (t)e-^ 



fceA ( ™> 



and the field W L °°\t,x) similarly defined, are space-homogeneous and partially 
(in the sense of the rotations of the torus) isotropic. 

Let H^' be the subspace of Hl correspondings to the modes with wavelengths 

in and consider the equation in H L n ' 

du^ + \vA L u { rt +n L n) B L (u^\u^A]dt = e ^ of d$ e~ ik ' x (B.3) 

(n) (n) 

where ir L ' is the orthogonal projection of Hl onto H L . 

Lemma B.l Ifu^ is a solution in Hl of AB.31 . with initial condition u^ n \0) and 
parameters (i/, L, 9), then 

u {n \t,x) := A^ (n) (A 1+/3 t,Ax) 
is a solution in Hl/\ of equation AB.31 with initial condition u^ n \0) and parame- 
ters ( z/A^ _1 , L/X, A z 13 9) (but with new Brownian motions). 



Proof. This statement is not clear a priori, especially because of the scaling 
transformation of the nonlinear term, so we give all the details. The solution u (n \ 
as a Fourier series, is given by 

u {n \t,x)= u {n \t,k)e- ik - x , 

and the solution u (n) , as a process in H^J^, is given by 

d n \t,x)= ^ u n \t,k)e- ik - x , 

The Fourier coefficients of u (n) and u (n) are related by the scaling 

u ( \t, k) = I u (n \t, x) e ik ' x dx = --r / u (n \X 1+l3 t, Ax) e ifc ' x dx 



Ld Jr L/x ' L d JT L/X 

(B.4) 

, k X '=J%) XP_ J u (n ){X l + ^ x , )e ik>. X < dx , = x ^ {X l + f3^ k , y 
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From the equation dB.3t in integral form, 

u (n) (t) + J \vA L vF> + tt^Bl (uW, uWJ] (s) ds 



ik-x 



fcGAV 



we have 



\P U ("\\1+P t , Xx) 



+ A 1+2 M vA L u^+^BM n \u^) (X 1+(3 s,\x)ds 



L/\aL/\, 



ik-x 



where j3^ X {t) := X + ^ /3^ x (X 1+f3 t) are new Brownian motions. The first term 

on the 1. h. s. is u^ n \t, x), and the first term on the r. h. s. is n (n ^(0, x). In addition, 
we have 

A L/X u^(t,x) = A 2+/3 (a l u^)(A 1+ ^,Ax). 
The proof of the claim will be complete if we show that 



A 



1+2/3 



A 1+ ^, Xx 



(t,x). 



For every ip G V L / X , by using the Fourier expression dB.21 ) of the non-linear term 
and the scaling of Fourier coefficients (IB.4I) . 



(^J x B L/x & n >M n >)(t, •), V )h l/x = (B L/X & n \ u (n >)(t, •), *T/x<P)B. 



» ~(n) 



(n) ~(n)\ 



» 



(l-u (t,h))u (t,l)-fi{k) 

h+l=k 

= X 1+2 ' 3 £ (W"W + ^)V n W + ^,{)- W 

= A 1+2 ^ L (u (m) , u<W)(A 1+ ^ A-), ^) Hl/x 
= X 1+2 P(^B L ( U (n) M n) )(X l+ % A-), ^ , 

where the sums above are extended to all wavelengths h, I and G ^-u x sucn that 
h + l = k. m 

B.2 Scaling theorem for stationary measures 

Similarly to section |2~.2.31 denote by V^ s (u, L, 9) the set of probability measures 
that are limit of homogeneous isotropic invariant measures of equations dB.3t . 



tin). 
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Given A > and f3 G R and /i G ^^5(^5 #)> let u be a random field on Tl 
with law define the random field u on T L /\ as 

n(x) = A^i^Aa;) 

and let Jl be the law of u on H L /x- More intrinsically, /I is defined by the relation 

/ f(u)djl(u)= [ /(A /3 n(A-))dMn) 
JH L/X Jh l 

for every bounded continuous / on H L j x . 



Theorem B.2 If y. G V% s (v,L,0) thenjle V% s , L/X, X~ 9 J. 

Proof. The measure \i of the theorem is the weak limit of a sequence {y nk } of 
invariant measures on Hj™ k ^ of the Galerkin problems with indexes rtfc. For each 
nfc, let u( n fc) be a stationary solution (on some probability space) of (IB.3I) . with 
parameters (y, L, 6) and marginal y nk . Let «( nfc ) be the rescaled process as above, 

which is a solution of dB.31 ) with parameters ^z/A^ _ 1 , Lj A, A (by the lemma 
above) and is a stationary process. Its marginal Jl nk is the scaling of y nk , similarly 
to the relation defined above between fi and Jl. Moreover Jl nk is an invariant mea- 
sure for equation dB.3b with parameters L/X, \~^~6^. From the weak 
convergence of \i nh to \i it is now easy to deduce the weak convergence of Jl nk to 
Jl. Therefore Jl G V^ s 1 , Lj A, A O^j . The proof is complete. ■ 
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